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Abstract 

We investigate the stability of the relativistic three-fermion system with a 
zero-range force in the light front form. In particular, introducing an in- 
variant cut-off, we study the dependence of the bound state on the coupling 
strength also for cases where the two-fermion system is unbound. The rel- 
ativistic Thomas collapse is discussed by solving the fully coupled integral 
equation system. Furthermore, we explicitly investigate the ground state 
mass of the three-fermion system and compare to previous simplified calcu- 
lations. 
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1. Introduction 

Relativistic constituent quark models in light front dynamics have re- 
ceived much attention lately [1]. A major advantage of this framework is 
the covariant description of few-body bound states under kinematical front 
form boosts [2]. This is a consequence of the stability of the light front 
Fock state decomposition under boost transformations [3]. Several studies 
have addressed nucleon properties, such as electromagnetic form factors, us- 
ing quarks as the relevant degrees of freedom and Gaussian wave functions 
on the light front [4, 5, 6]. Other calculations of the proton electric form 
factors [1] use nucleon null-plane wave functions obtained from the solution 
of the Faddeev equation [7] with a zero-range force [8] acting between the 
constituent quarks. In this effective model the spin degrees of freedom and 
confinement are missing. Naturally, zero-range interactions provide a simple, 
but important limiting case for short range forces and there has been renewed 
interest in zero-range models with applications in nuclear and atomic physics 
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in the past years [9, 10, 11, 12]. The present work is an extension of the devel- 
opments given in Ref. [8], where the three-boson bound state problem on the 
light front was considered. Subsequently, the relativistic three-boson problem 
with zero-range interactions using light front dynamics has been reformulated 
by Carbonell and Karmanov in a covariant way [11]. While in the former 
publication the zero-range interaction was smeared, the authors of [11] only 
assumed the existence of a two-body bound state for the UV regular izat ion 
without introducing further restrictions in the Faddeev equation itself. With 
the latter regularization scheme, which is called M^-regularization in the fol- 
lowing, the mass of the three-body bound state may vanish despite the finite 
value of the mass of the bounded two-body system. This is the relativistic 
analog of the so-called Thomas collapse. Originally, it was pointed out that 
non-relativistic three-body systems, based on zero-range forces, experience 
a Thomas collapse [13]. Since the inverse of the binding energy | J5 3 1 ~ 1 is a 
measure of the size of the system the three-body state collapses when the 
binding energy is unbounded from below, i.e. B 3 — > — oo. 

A flaw of the M 2 £-regularization scheme, as used in [11], is that the ex- 
plicit dependence of the three-body mass on the coupling constant is not 
recovered, since the three-body bound state mass is parameterized by the 
two-body mass. In order to investigate the three-body bound state equation 
even if no two-body bound state exists, the introduction of a different regu- 
larization scheme is necessary. To this end an invariant cut-off A as proposed 
by Lepage and Brodsky in [14] was introduced into the three-boson equation 
in Ref. [15] . More specifically, the requirement is that the masses of the two- 
and three-body systems are smaller than the cut-off, i.e. M| , M| < A 2 . It 
has been shown in Ref. [15] that in the Lepage-Brodsky (LB) regularization 
scheme the limit of very large A coincides with the results of Carbonell and 
Karmanov where a comparison can be made. In particular, the relativistic 
Thomas collapse appears for any value of A. 

The validity of the bosonic model is evident by the calculation of the 
proton electric form factor obtained from the Faddeev wave functions, that 
reproduces well the experimental data for low momentum [1]. The dominance 
of this scalar channel for the description of the nucleon static properties has 
been confirmed by systematic investigations of the relativistic quark coupling 
effects in the nucleon electromagnetic form factors [16]. The scalar coupling 
between the quark-pair is preferred by neutron data, independently on the 
choice of the momentum wave functions [17, 18, 19]. 

However, to achieve a more realistic description of the nucleon the spin 
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degrees of freedom have to be accounted for. This leads to a coupled system 
of integral equations instead of one Faddeev equation. The first step in this 
direction has been done again by Karmanov and Carbonell [20], where the 
first results using the M2B- regularization for a two-fermion kernel have been 
reported. The well-known difficulties to treat the spin on the light front can 
be reduced by using kernels factorized relative to initial and final states. In 
particular, Karmanov and Carbonell have chosen a simplified kernel with 
properties similar to the bosonic interaction. Indeed, the successes of the 
bosonic model describing the proton properties indicate that this limitation 
is reasonable. 

In this paper we revisit the interaction given in [20] using the LB regu- 
larization scheme. In doing so we relax their limitations and investigate 

(i) the three-fermion equation even if no two-fermion subsystem exist, 

(%%) the direct dependence of the two- and three-body bound states on the 
coupling strength, 

(Hi) in contrast to [20], the full set of coupled integral equations and 

(iv) the critical coupling for the Thomas collapse. 

This paper is organized as follows. In the Section 2 we present the deriva- 
tion for the two- and three-fermion bound state equation on the light front 
using a general zero-range interaction where the kernel is factorized relative 
to the initial and final states. We consider a specific interaction kernel cor- 
responding to a 1 So nucleon in the Section 3 and derive the coupled integral 
Faddeev equations. Our results are compared to the previous calculations of 
Carbonell and Karmanov and to the three-boson case in the Section 4. A 
discussion of the main findings and perspectives is given in the conclusions. 

2. 2-Fermion and 3-Fermion bound states for separable kernels 

The zero-range two-fermion kernel can be constructed utilizing different 
spin couplings. We use kernels being factorized relative to initial and final 
states, thus having the form 

/C&| (1, 2; 1', 2') = XK aia2 (1, 2) K« (l', 2') , (1) 

with 

K aia2 (1, 2) oc [u ai (px) ufiVu n (p 2 )] (2) 
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Figure 1: Equation for the two-body i-matrix with zero-range interaction. 

where V indicates the interaction vertex, <7j denote the spin indices of the 
particles and u ai the corresponding spinors. The factorization (of the kernel) 
leads to simplifications that are independent from the specific form of the 
spin coupling. 

2.1. The two -body problem 

We briefly review the basic ingredients of the two-body equations on the 
light front given previously in [1, 8, 15]. For a zero-range interaction the 
equation represented by Fig. 1 can be summed leading to a solution for the 
two-particle propagator t(P 2 ), i.e. 



t(P 2 ) = (iX- 1 - B(P 2 ))~ 1 , 



(3) 



where P 2 is energy-momentum vector of the two-body system. The expres- 
sion for B(P 2 ) corresponds to a loop diagram and is given by 



B{P2) 



dx'd 2 k' 



K 



x'{l-x')Pi-M^ 



(4) 



where 



K = 



(k'l + m 2 )/x'(l-x') 
A^Trj/C}. 



(5) 
(6) 



M 2 q is the mass of the virtual two-particle state and m the mass of the 
constituent particles. The integration variables in the loop are x', k' ± are 
the relative coordinate of the two-body system as defined in Appendix A. 
Equation (4) is a generalization of the corresponding equation in the bosonic 
case [1, 8, 15], where the spin structure has been missing. This which is 
equivalent to K — 1 in eq. (4). The integral (4) has for a contact inter- 
action an UV divergence that has to be regularized. If the kernel leads to 
a logarithmic divergence, like for the scalar bosons, it can be absorbed in a 
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redefinition of A by assuming that the two-particle propagator t(P 2 ) has a 
pole at P 2 = M| B . Then holds the relation 

iX- 1 = B(M 2B ). (7) 

Hence, the renormalization scale of the propagator is fixed by the mass of the 
two-particle M 2 b- With the M 2 £-regularization the condition (7) imposed in 
the denominator of eq. (3), i.e. B(M 2 b) — B(P 2 ), cancels the UV divergency 
of t(P 2 ). The assumption of the existence of a two-body bound state reg- 
ularizes the loop integral completely only for a logarithmic divergence. For 
completeness, we remark that, in general, additional conditions have to be 
introduced in order to fully regularize the self-energy B(P 2 ). Note, that us- 
ing the M 2 £-regularization scheme the relation between the two-body mass 
and the coupling constant is not resolved. The two-particle propagator is the 
input for the relativistic three-body equations. In order to investigate the 
three-body bound state equation also for cases where no two-body bound 
state exists another regularization scheme has to be used. This allows us to 
find the relation between the coupling constant and the two- and three-body 
masses. 

One can introduce an invariant cut-off A in the integral (4). The require- 
ment that the invariant mass of the two-particle system is smaller than A 
(LB regularization), i.e. 

Mf < A 2 , (8) 
makes the integral (4) finite. Thus, the integral limits are 



J dx' j d 2 k' ± 



^2tt / / k' ± dk' ± , (9) 



x' • 

mm 



where 



x 

max 



X - (l - y/1 - 4m 2 /A 2 ) , (10) 

\ (l + y/1 - Amy A") , (11) 
KL X = AV(l-x')-m 2 . (12) 

Consequently, t(P 2 ) — > t\(P 2 ) depends on A as well. The LB regulariza- 
tion scheme has been utilized in [15] and the subsequent calculations in [21]. 
Accordingly, one can calculate the pole of the two-particle-propagator deter- 
mining the bound state mass M 2B for any value of the cut-off A. 
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Figure 2: Diagrammatic representation of the Faddeev equation for a zero-range interac- 
tion, eq. (15). The gray area indicates the two-body input given in Fig. 1 and eq. (3). 



2.2. The three-body problem 

Investigating the three-body system, we consider the three-body bound 
state wave function ^(1,2,3) on the light front, where 1,2,3 shorthand for 
Pi,P2,P3- The wave function is related to the vertex function T by 

tt(l,2,3) = m*-mI b ^ - M2 = (^+^+P3) 2 , (13) 

where M 3 b is the three-body bound state mass. We expand the vertex func- 
tion T in the standard way as 

r(i, 2, 3) = r 12 (i, 2, 3) + r 23 (i, 2, 3) + r 3 i(i, 2, 3), (14) 

where the Faddeev components 1^(1,2,3) satisfy equations involving the 
interaction between particles (ij) only. For three fermions the equation for 
Faddeev component r 12 (l,2,3), by adding the spin indices, reads 

d 2 k' L dx' 1 



C(l, 2 ,3)^/ ^/C^(l,2;l',20^ 



- x>) M| - Pi 



x 



K,,(l', 2', 3) + T°,, (2', 3, 1') + K,,, (3, 1', 2' 



(15) 



where a is the spin projection of the bound state. A diagrammatic rep- 
resentation of this equation is given in Fig. 2. In general, as explained in 
Appendix A, the arguments of r^ 10 . (1, 2, 3) can be rewritten in variables 
kij_ and Xi (i = 1, 2, 3), where the transversal momenta are ki± = pi± — XiP^± 
and the fraction of the total plus component carried by the single particles 
is given by X{ = pf / P% ■ Here P^ denotes the momentum of the three-body 
system. Alternatively, one can rearrange the now defined variables of two 
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particles % and j to the corresponding pair variables kj_ and x introduced 
for the two-body problem. This is very useful for eq. (15), where a natural 
split to a pair (ij) and to a spectator particle (the third one) emerges. The 
expression for the two-body mass embedded in the three-body system is 

Pi = (Ps-Psf 

= (M 3 -pi)(^M 3 - P ^±^-p 3 l, (16) 

where the momentum of the three-body system is taken at rest, i.e. P 3 = 
(M 3 , M 3 , 0,0) in standard light front coordinates. Using the variables intro- 
duced the following relation 

p2 /, Uyf2 h 2 ± + (1 - x 3 )m 2 

P 2 = (1 - x 3 )M 3 , (17) 

x 3 

holds. 

In this work we investigate a bound state with total spin a = 1/2. Because 
of the factorized form of the kernels, the Faddeev component r£ 1(7 (1, 2, 3) 
must be proportional to K aia2 (1, 2) and hence has the form 

(1, 2, 3) = K aia2 (1, 2) Gl z (1, 2, 3). (18) 

Using eq.(15) one finds the equation for G° (I, 2, 3) 

«:,(i.2.3) = ^ F / E^ 1 - 2 ; 1 ' 2 W-*')M1,-P| 

x (1', 2')G- (!', 2', 3) + (2', 3)G% (2', 3,1') (19) 

+ ^ i (3,l / )^,(3,l',2 / ) 

The functions G" are the components of a 2 x 2 matrix labeled by the spin in- 
dices. In order to decompose G° 3 into basis structures we follow Refs. [22, 23], 
i.e. represent it as a projection to spin states of four independent elements 
from the Dirac space. Introducing the orthogonal 4x4 matrices {Si} the 
matrix G a az can be written as 

4 

(1,2,3) = Y J 9i{^^)uMS i u"{P 3 ). (20) 
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The explicit formulae for the {Si} are given in the Appendix B. Note that 
the basis elements satisfy the following orthogonality conditions 



Tr 



with Si = 70^70. The matrices are not normalized to one rather 



Sty +m)Sjtf 3 + M 3 







for i ^ j, 



(21) 



Tr 



Sity 3 + m)Sitf 3 + M 3 ) 



Ni 



8x 3 k 2 
8x 3 m 2 
8x 3 k 2 C 2 s 
8x 3 m 2 C 2 



for % — 1, 
for i = 2, 
for i — 3, 
for i — 4. 



(22) 



The factors C ps and Cp S are given by [24] 
1 



m 2 u; ■ P 3 



C l = — — 

ps m 4 



(23) 

where u is a four-vector with u 2 = determining the light front orienta- 
tion by the equation uo ■ x — 0. We use w = (1,0,0,-1) to recover the 
standard LF approach. We keep u unspecified in the formulae (23) for short- 
ness. Observe that the matrix K aia2 (l,2) in eq. (18) is antisymmetric with 
respect to the permutation 1 f> 2. The functions #1,2 are symmetric, i.e. 
#1,2(1,2,3) = #1,2(2, 1,3). On the other hand the function #3,4 are antisym- 
metric, i.e. #34(1,2,3) = —#34(2,1,3), because of the additional antisym- 
metric factor C ps . Substitution of the decomposition (20) into eq. (19) leads 
to 



^#,(1,2,3) 



A 



d /c 1 doc 



1 



(2tt) 3 J 2x'(l - x') Mf - P 



(24) 



x 



4 4 
KN' i9 i{l', 2', 3) + <^( 2 ', 3, 1') + v ii9i&, ^ 2 



The factors N' 12 are the same as Aq )2 , whereas in JVg 4 the factor Cp S is 
replaced by 



m 2 w • P3 



e^Pi,P2 v P 3p ^e^p'p^P 3 ^ 



((k' ± ■ k±)kl ± - (k± ■ k 3 ±)(k' ± ■ k 3± )} 

TYli V / 



(25) 



Furthermore the kernels and i>?- are given by 

<• = y2u„(P 3 )S t u^(p 3 )K^ (l',2')K a , a3 (2',3)u a[ (p' 1 )S' b ^(P 3 ), 

±( ( 26 ) 
4 = ^u^S^fo)^ (l',2')^ (3,l>^(p 2 )^(P 3 ). 

In Appendix B we explicitly explain how the components of S' b and S' c are 
constructed. The factorization of the contact kernel ensures that g% = = 0, 
as shown in Ref. [22]. Similarly one proofs that gi and g 2 depend on k 3 ± and 
x% only, which we rename to q± and y respectively. Therefore the contribution 
of v\j and i>?- are equal to each other. Additionally, since the first term g i 
on the r.h.s. of eq. (24) does not depend on the integration variables, we 
can rearrange the integral equation in terms of the two-particle propagator 
t(P 2 ) given in eqs.(3)-(4). The general form of the integral equations for the 
three-body bound state reads 

7Vi(27r) 3 7o J ^ ± + x'q ± ) 2 + m 2 -x'(l-x')P 2 

= j(p^_ f dx , r^^yd-^^D+^^d-.)^!) 

N 2 (2nf J J L (k' ± + x'q ± ) 2 + m 2 -x'(l-x')P 2 

(27) 

Finally, replacing x'(l — y) — > x' the equations (27) is transformed to 
t(P 2 ) f'- y J , f j27 , ViMx',k' ± )+v 12 g 2 (x',k' ± ) 

(28) 



t(P 2 ) f y f 
9l{y ^ ) = N^fj dX 'J d2k 

a (v a, ) - m) r V dx> fd 2 k> ^9,{x' A) + v 22 g 2 {x' A) 



where Mg is the mass square of the virtual three-particle state in the rest 
system, i.e. 

m2 = k'j + m 2 + q 2 + m 2 + (g + q)\ + m 2 
30 x' y 1 — x' — y 

In order to consistently perform the LB regularization in the two- as well in 
the three-fermion equations we also constrain the mass square of the virtual 
three-particle state by 

J " dx' j d 2 k' ± ->• J V dx' j d 2 k' ± 6(M 2 - A 2 ). (30) 



After fixing an explicit interaction kernel and spin structure the aim is 
to solve the eigenvalue problem posed by eqs. (28) in order to obtain the 
mass Mf = P 3 2 of the three-fermion system as a function of the the coupling 
constant A. 



3. The 1 S Model 



In this section we investigate a particular kernel which leads to significant 
simplifications in the three-body equations (28). The spin structure of this 
kernel is given by 



K, 



<Tl<72 



;i,2) 



m 



2') 



w • (pi + p 2 ) 
m 

u ■ (pi + p 2 ) 



Ki (Pi) i^U c u a2 (p 2 )] , 



^ (p 2 ) C/ c ^ 7 5^ (Pi) 



(31) 



where C/ c = ?7°7 2 denotes the charge conjugation operator. In the nonre- 
latististic limit the kernel reads 



«'(l,2;l',2') oc ^(<mXJK*w) 



(32) 



= AC* 



oo 

(l/2)(7i(l/2)<r 2 



C(T/2)ai (1/2)4 ("'^^O (^4 i(7 2/^i) > 



where w a represent the two-component (Pauli) spinors and CR 



(1/2)cti(1/2)ct 2 



are 



the Clebsh-Gordon coefficients [25]. Therefore, this kernel corresponds to an 
interaction in the 1 S state. 
For this interaction we find 



K = \- lr Tr{)C} 



m 



J2 (^2) U c t0 l5 u< (pi)] (pi) ^75^ (p' 2 )] 



(w • P 2 ) 2 ^ 
= 8\m 2 x(l — x). 
From equation (4) the loop integral S 2 (-P 2 ) reads 

4m 2 



B(P 2 ) = -i- 



(2tt) 3 

Using the M 2 b regularization one finds 



/ dx'cfk'^-^- 
J * 2 



2 -^2 2 ' 



48tt 2 
10 



(33) 



(34) 



(35) 



where 



k(P 2 , M 2B ) 



3T 2 (P> 



with 



T 3 (P 2 ) 
3T 2 (M 2B ) + 1 
T 3 (M 2B ) 

_ P2 



arctan Y(P 2 



T 2 (P 2 ) 



arctan T(M 2B ) + 



1 



T(P 2 ) = 



y/4m 2 -P; 



T 2 (M 2B ) 

for < P| < 4m 2 , 
for P| < 



(36) 



(37) 



v /4m 2 -P| 

Introducing the invariant cut-off as in eq.(8) the loop integral can be 
directly evaluated using the limits of the integration given in eqs.(10)-(12). 
In this way, for any value of the cut-off and the coupling we can calculate 
the pole of the two-fermion propagator. 

Now we turn to the three-fermion case. Calculating the 2x2 matrix for the 
kernel (31) using eq.(26) yields 

v u = 32m 2 yx' 2 (f ± + k' ± - q ± ^j , 

V12 = 0, 

i>2i = 32m 4 yx' (1 — 2x') 
v 11 = — 32m 4 yx' (1 — x') . 

Shifting x'{l — y) — > x' and using z = x'/ (1 — y) we find the following set of 
equations 

4m 2 t(P 2 ) f l ' y dx' 



(38) 



9i(y,q±) = 



92(y,q±) = 



(27T) 

(l + k' ± - q ± /g 2 )z 2 

n- 

AmH(P 2 ) 



x'(l — y — x') 



d 2 k' ± 9(Mi - A 2 ) 



1-2/ 



(27T) 



9i (x',k' ± ), 
dx' 

-x') 



(39) 



x'(l — y 



d 2 k'J{M 2 - A 2 ) 



z(l - 2z) 9l (x', k' ± ) - z(l - z)g 2 {x', k' ± ) 



p2 



Ml 



(40) 



The scalar function g 2 enters in the second equation only. The first equa- 
tion differs from the three-boson equation by the factor (1 + k' ± • q±/q±)z 2 . 
We point out that the complete system of equations must be solved, not only 
the first equation as in Ref. [20]. 
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Figure 3: The solution of the two- (dashed) and three- fermion (solid) bound state equations 
as a function of the effective coupling A e g = Am 2 . Horizontal dashed line shows the two- 
fermion break-up. 

4. Results 

The solution of the equations (39)- (40) allows us to determine the three- 
fermion mass. In Ref. [20] eq. (39) is solved using the M 2 b regularization. 
Instead we consider the LB regularized version of (39)- (40). In order to 
compare our results with those of [20] we set the cut-off to a large but finite 
value, i.e. A = 10 15 m, see Ref. [15]. For the bosonic system it was shown 
by one of the authors that for large invariant cut-off the resulting function 
M 3 b(M 2 b) coincides with the analogous results in M 2 b regularization [15, 21]. 
Therefore, we expect this feature as well in the fermionic case. 

Fig. 3 shows the masses M 2 b of the two-fermion bound state and M 3 b of 
the three-fermion bound state as a function of the effective coupling strength 
A e ff = Am 2 . All masses M 2 b and M%b are given in units of the elementary 
fermion mass m. 

Note the following features: 

(i) The two- as well as the three-fermion mass are monotonic function of 
the coupling. In particular, for the critical coupling X c eS w 6.972 the 
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three-fermion bound state mass vanishes. This is the relativistic analog 
of the Thomas collapse [13] and is known to occur for the three-boson 
bound state with a similar zero-range interaction, see e.g. [15, 21]. 
This critical coupling corresponds to a critical value of two-fermion 
mass M£g ps 1.42m. In comparison, the critical two-fermion mass 
M| B 1.35m calculated in [20] is notable lower. 

(ii) As M 2B 2m, the three-fermion bound state still exists. This may 
give rise to the Efimov effect [26, 27], an interesting issue for further 
investigations. 

(Hi) There is a region of parameter space where both M 2 b and M 3 b exist. 
Only in this case it is possible to plot M 3B against M 2B . 

We focus on the investigation in the region of the coupling strength A, 
where the two-fermion bound state exists. Therefore, by comparison to the 
results of [20], we can estimate the role of the coupled structure of (39)- (40). 
In Fig. 4 the mass M 3B is plotted against M 2B . We have also included the 
solid line for fermions and the dash-dotted line for bosons. For comparison 
to the calculation of [20] we independently solve, using the LB regularization 
scheme, eq. (39) solely and present the results in Fig. 4 as dashed line. 
Furthermore, the dotted line indicates the threshold for the three-body mass. 

The resulting functions M 3B (M 2B ) for the bosonic and for the fermionic 
cases have the same qualitative behavior. By decreasing M 2B the three-body 
mass decreases very quickly and vanishes at a critical values M 2B of the 
two body mass. The value of the critical mass for the bosonic system is 
M 2 b ~ l-45m and greater than the corresponding values for the fermionic 
cases, both for the coupled equation system and for the solution of the first 
equation only. 

It is interesting to see why there is this notable discrepancy between the 
whole calculation and the result obtained by solving (39) alone. 

The latter method automatically excludes the trivial solution g 1 = for 
eq. (39). However, this solution is physically non-trivial in the coupled 
system. 

In more detail, the system of integral equations (39)- (40) can be written 



as 




(41) 
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Figure 4: Three-fermion bound state as a function of M 2 b (solid line) in comparison to the 
three-boson bound state (dash-dotted line) and to three-fermion bound state calculated 
by solving equation (39) alone, with A —> oo (dashed line). To completeness, the threshold 
for the three-body mass is indicated (dotted line). 

where the operator V collects all contributions from and I is the identity 
operator. 

To achieve a numerical solution we have discretized the equations (39)- 
(40) using a momentum space grid 7&, with k = 1, ■ • • ,n, where n is the 
number of the lattice points. Each of the scalar function gi, % = 1,2 be- 
comes a vector g\ of dimension n. The operator eigenvalue problem (41) is 
approximated by the corresponding matrix problem 



All possible solutions fulfill of course the condition det(VK) = 0. For the 1 S'o 
model (31) one has W12 = and the solvability condition reads 



W g = with 




(42) 



where 




(43) 




(44) 
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Figure 5: Mass of the three-fermion ground state (solid line) and of the first three- fermion 
excited state (dashed line) calculated by solving the equation system (39)- (40) as a function 
of the effective strength A g . 

Evidently, a solution of the whole system satisfies det(Wn) = or det(W22) = 
0, whereas the solutions of eq.(39) satisfies det(Wn) = only. Hence, solving 
solely eq.(39) the physical solutions corresponding to the case g\ = can be 
lost. 

Because for arbitrary coupling strength A there exist more solutions to 
(42), we consider the two smallest eigenvalues. In Fig. 5 the mass of the 
ground state (lowest eigenstate) and the first excited state (next-to-lowest 
eigenstate) depending on effective coupling strength X eS are plotted as solid 
and dashed line, respectively. The ground state coincides with the solution 
shown in Fig. 3. For the first excited state we find that the Thomas col- 
lapse occurs at X e g ~ 6.98 that corresponds to a critical two-fermion mass 
M|g « 1.35m which is in agreement with the critical value obtained by solv- 
ing eq.(39) only We therefore suspect that the solution given in Ref. [20] 
describes the first excited state. This observation is confirmed in Fig. 6 
where the bound state mass solving only eq. (39) (crosses) is compared with 
the solution for the first excited state (circles) obtained by solving the cou- 
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X 



eff 



Figure 6: Mass of the three- fermion bound state calculated by solving eq.(39) only (crosses) 
and mass of the first three-fermion excited state calculated by solving the set of equa- 
tions (39)-(40) (circles) as a function of the effective strength A e ff . 

pled equation system (39)- (40) for different values of the effective coupling 
strength A e fj. The two masses agree well for all A e fj but small derivations 
are noticeable close to the breakup, i.e. for A e ff < 6.86. In other words 
the solution of eq.(39) describes three-fermion bound states, but the ground 
state is missing. The true ground state is only obtained in the solution of 
the coupled system of integral equations. 

5. Conclusions 

We have derived and solved the coupled equations of a relativistic three 
fermion system subject to an effective scalar zero-rage interaction in the 1 S'o 
channel. In the computation of the three- and two-fermion bound states as 
function of the coupling constant A we have introduced an invariant cut-off 
A (LB regularization scheme) and solved the equations for large A. Addi- 
tionally, we have investigated the solution of eq. (39) alone with the same 
regularization scheme. Thereby the results of Ref. [20] were confirmed. Anal- 
ogous to the three-boson system we find the relativistic Thomas collapse 
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in both calculations. In particular, solving the coupled equations we find 
that the two-fermion critical mass M| B ps 1.42m is different from the value 
given in [20]. This is explained by the fact that in the computation of Car- 
bonell/Karmanov the lowest three- fermion mass is missing. Instead of the 
ground state their computed lowest mass seems to agree with the mass of 
the first excited state. Taking in account the full equation system, we have 
obtained a critical two-body mass different from the one in the literature and 
more close to the critical value in the bosonic case (see Fig. 4). 

In particular, our results provide evidence we that the bosonic or fermionic 
nature of the constituents does not play a very important role in this model. 
Consequently, using the bosonic-type Faddeev equations to describe the iso- 
lated nucleon [15] or embedded in a quark matter medium [28, 29] can be 
justified a posteriori. 

Nevertheless, an implementation of the spin degrees of freedom is de- 
sirable and a first consistent step in this direction is studied here. Further 
investigations of the three-fermion correlations in hot and dense quark matter 
following [29] can be performed using the interaction kernel (31). 

The theoretical methods presented here are model independent and allow 
the investigation of the nucleon with more complicated kernels, such as the 
one corresponding to an interaction in the P-wave state with S = 1 [22] . 

Acknowledgment: The authors are grateful to V.A. Karmanov for the 
providing his theoretical notes [22]. This work is supported by Deutsche 
Forschungsgemeinschaft and by the Helmholtz International Center for FAIR 
within the LOEWE program of the State of Hesse. 

A. Kinematics 

Let a M = (a , a 1 , a 2 , a 3 ) be a generic four vector. The light front dynamics 
is given by the transformation 

a?" = (a + , a~ , a 1 , a 2 ) = (a + ,a~,a_i_) , (45) 

where 

a ± = a°± a 3 . (46) 

The covariant components follow from a M = g^ v a v , where the metric 
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tensor is given as, 
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(47) 



The mass shell condition for a particle of mass m and momentum p reads 

(48) 



p + p 



-2 i 2 

p ± + m . 



For the description of the kinematics of two- and three-particle systems 
it is convenient to introduce the total momentum of the particles and some 
internal variables. For two particles with momenta p\ and p2 and masses mi 
and m 2 respectively we define 



Pi = P1+P2, 
x = pf/P+, l-x = pt/P 2 + , 
k = (1 - x)pi - xp 2 . 



(49) 
(50) 
(51) 



The invariant energy square, Si 2 , i.e. the mass of the virtual two-body state, 
is given by 



S12 = Mf 



p+ p- _ p2 
r 2 r 2 r 2± 



2 H2 2 
Pi± + m i 



E 

i=i 



kt 



2_L 



x 1 — X 



2± 



1 — rr x ^ x(l — x) 



(52) 

(53) 
(54) 



Choosing P 2 ± = and for equal masses, i.e. m\ = m 2 = m we recover the 
expression in eq.(5). 

For three particles with momenta pi, p 2 and ^3 and masses mi, m 2 and 
m 3 we analogously introduce the plus momentum fractions and the total 
momentum as 



Xi 



Pi +P2+P3, 

3 

P+/P+, 5> = 1. 



(55) 
(56) 
(57) 
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In general one can introduce Jacobi coordinates to treat the kinematics of the 
three-body system [30]. However, because of the structure of the Faddeev- 
equation, where a spectator particle and a pair of two particles emerge, we 
prefer to express the kinematics using the variable of the spectator particle 
-say the third particle - and the one of the pair (12), where kinematics has 
been discussed above. Note that the momentum fraction of the three-body 
system X; L and the variable x do not coincide but are related by 

xi = x(l-x 3 ), (58) 
x 2 = (l-x)(l-x 3 ). (59) 

The additional constraint used in this paper, namely that the three-body 
system is at rest, allows us to connect the single transverse momenta of the 
particle in the pair with the relative transverse momentum k± by 

Pi± = ~k± - xp 3± , (60) 
P2± = k ± - (1 -x)p 3J _. (61) 

The invariant energy square, s 12 3, i.e. the mass of the virtual three-body 
state, can be written as 

Sl23 = M 3 2 = P 3 + P 3 -PL = i:^T^ P s- p L (62) 



i k i 



3 H2 2 

Pi± + < 52 



i=l X i 



~ PL- (63) 



Demanding P 3 ± = and for equal masses, i.e. mi = vri2 = m 3 = m we 
recover, with the necessary variable identification, eq.(29). 

B. Orthogonal basis 

In this section we define the orthogonal basis {Si} as given in Refs. [22, 23] 
used for the decomposition of the 2 x 2-matrix G° 3 in eq. (20). We choose 
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the form 



Si 



2x 3 - (m + x 3 M 3 ) 



u-P 3 



m 







2x 3 - (m - x 3 M 3 







w-Ps 



«Cps75, 



to- 







u-P 3 

where the factor C ps is given in eq.(23). 
Furthermore, the matrices Si = 7o<S'i7o explicitly read 



(64) 
(65) 
(66) 
(67) 



5i 

S 3 

s 4 



2x 3 - (m + x 3 M 3 ) 







w-P 3 







2rr 3 - (to - X3M3) 







w-P 3 







(68) 
(69) 
(70) 
(71) 



Note that 5i )2 = Si )2 . These matrices satisfy the orthogonality (21) and 
normalization conditions (22). 

The permutated structures S' b and S' c are derived from eq. (64-67). S' b 
is constructed using the permutation (1' — > 2' — > 3 — > 1'), e.g. for S' bl one 
performs the substitution x 3 — >■ x[ = (1 — x)(l — x 3 ) in eq. (64), while for 
S b3 and S' M one additionally has to replace C ps by 

1 „,™ , , „ 1 



CL = 



(72) 



in eq. (66) and (67) respectively. Analogously, S' c is derived from S by the 
permutation (1' — > 3 — > 2' — > 1'), e.g. for S^i the substitution x 3 — >■ x 2 = 
x(l — £3) has to be used in eq. (64), while for S' c3 and S^ 4 one has to use 

C - = SS^""^^ - ^^V-V, (73) 
instead of C ps in eq. (66) and (67) respectively. 



20 



References 

[1] W. R. B. de Araujo, J. P. B. C. de Melo, T. Frederico, Faddeev null 
plane model of the nucleon, Phys. Rev. C52 (1995) 2733-2737, doi: 
10.1103/PhysRevC.52.2733. 

[2] B. D. Keister, W. N. Polyzou, Relativistic Hamiltonian dynamics in 
nuclear and particle physics, Adv. Nucl. Phys. 20 (1991) 225-479. 

[3] R. J. Perry, A. Harindranath, K. G. Wilson, Light front Tamm- 
Dancoff field theory, Phys. Rev. Lett. 65 (1990) 2959-2962, doi: 
10.1103/PhysRevLett.65.2959. 

[4] I. G. Aznaurian, A. S. Bagdasaryan, N. L. Ter-Isaakian, Relativistic 
quark model in the infinite momentum frame and static characteris- 
tics of nucleons, Phys. Lett. B112 (1982) 393-396, doi:10. 1016/0370- 
2693(82)91076-0. 

[5] Z. Dziembowski, Relativistic model of nucleon; static properties and 
electromagnetic soft form-factors, Phys. Rev. D37 (1988) 778, doi: 
10.1103/PhysRevD.37.778. 

[6] P. L. Chung, F. Coester, Relativistic constituent quark model 
of nucleon form- factors, Phys. Rev. D44 (1991) 229-241, doi: 
10.1103/PhysRevD.44.229. 

[7] L. D. Faddeev, Scattering theory for a three particle system, Sov. Phys. 
JETP 12 (1961) 1014-1019. 

[8] T. Frederico, Null plane model of three bosons with zero range interac- 
tion, Phys. Lett. B282 (1992) 409-414, doi:10. 1016/0370-2693(92)90661- 
M. 

[9] D. V. Fedorov, A. S. Jensen, Regularized zero-range model and an ap- 
plication to the triton and the hypertriton, Nucl. Phys. A697 (2002) 
783-801, doi:10.1016/S0375-9474(01)01266-0. 

[10] D. V. Fedorov, A. S. Jensen, Correlation induced collapse of many-body 
systems with zero-range potentials, Phys. Rev. A63 (2001) 063608, doi: 
10.1103/PhysRevA.63.063608. 



21 



[11] J. Carbonell, V. A. Karmanov, Three-boson relativistic bound states 
with zero-range interaction, Phys. Rev. C67 (2003) 037001, doi: 
10.1103/PhysRevC.67.037001. 

[12] M. T. Yamashita, T. Frederico, L. Tomio, Neutron— 19 C scatter- 
ing near an Efimov state, Phys. Lett. B670 (2008) 49-54, doi: 
10.1016/j.physletb.2008.10.030. 

[13] L. H. Thomas, The Interaction Between a Neutron and a Proton 
and the Structure of H**3, Phys. Rev. 47 (1935) 903-909, doi: 
10.1103/PhysRev.47.903. 

[14] G. P. Lepage, S. J. Brodsky, Exclusive Processes in Perturba- 
tive Quantum Chromodynamics, Phys. Rev. D22 (1980) 2157, doi: 
10.1103/PhysRevD.22.2157. 

[15] M. Beyer, S. Mattiello, T. Frederico, H. J. Weber, On the stability of 
three-body bound states on the light front, Few Body Syst. 33 (2003) 
89-97, doi:10.1007/s00601-003-0014-2. 

[16] W. R. B. de Araujo, E. F. Suisso, T. Frederico, M. Beyer, H. J. We- 
ber, Relativistic quark spin coupling effects in the nucleon electromag- 
netic form factors, Phys. Lett. B478 (2000) 86-93, doi:10.1016/S0370- 
2693(00)00238-0. 

[17] W. R. B. de Araujo, T. Frederico, M. Beyer, H. J. Weber, Neutron charge 
radius: Relativistic effects and the Foldy term, Int. J. Mod. Phys. A18 
(2003) 5767-5780, doi:10.1142/S0217751X03015866. 

[18] E. F. Suisso, W. R. B. de Araujo, T. Frederico, M. Beyer, H. J. We- 
ber, Relativistic quark spin coupling effects in the correlations between 
nucleon electroweak properties, Nucl. Phys. A694 (2001) 351-371, doi: 
10.1016/S0375-9474(01)00978-2. 

[19] W. R. B. de Araujo, T. Frederico, M. Beyer, H. J. Weber, mu(p) 
G(E(p))/G(M(p)) and qF2(p)/Fl(p) in a relativistic quark model, Eur. 
Phys. J. A29 (2006) 227-234, doi:10.1140/epja/i2006-10072-2. 

[20] V. A. Karmanov, J. Carbonell, Critical stability of three-body relativis- 
tic bound states with zero-range interaction, Few Body Syst. 34 (2004) 
85-90. 



22 



[21] S. Mattiello, The stability of the relativistic three-body system and in- 
medium equations, Few Body Syst. 34 (2004) 119-125. 



[22] V. A. Karmanov, Three-fermion relativistic bound states with zero- 
range interaction, private communication (2004) . 

[23] V. A. Karmanov, J. F. Mathiot, A. V. Smirnov, Nonperturbative cal- 
culation of the anomalous magnetic moment in the Yukawa model 
within truncated Fock space, Phys. Rev. D82 (2010) 056010, doi: 
10.1103/PhysRevD.82.056010. 

[24] V. A. Karmanov, The nucleon wave function in light-front dynamics, 
Nucl. Phys. A644 (1998) 165-200, doi:10.1016/S0375-9474(98)00579-X. 

[25] J. Carbonell, B. Desplanques, V. A. Karmanov, J. F. Mathiot, Explicitly 
covariant light-front dynamics and relativistic few-body systems, Phys. 
Rept. 300 (1998) 215-347, doi:10.1016/S0370-1573(97)00090-2. 

[26] V. Efimov, Energy levels arising form the resonant two-body forces in a 
three-body system, Phys. Lett. B33 (1970) 563-564. 

[27] E. Braaten, H. W. Hammer, Universality in Few-body Systems 
with Large Scattering Length, Phys. Rept. 428 (2006) 259-390, doi: 
10.1016/j.physrep.2006.03.001. 

[28] M. Beyer, S. Mattiello, T. Frederico, H. J. Weber, Three-quark clusters 
at finite temperatures and densities, Phys. Lett. B521 (2001) 33-41, 
doi:10.1016/S0370-2693(01)01175-3. 

[29] S. Strauss, S. Mattiello, M. Beyer, Light-front Nambu-Jona-Lasinio 
model at finite temperature and density, J. Phys. G36 (2009) 085006, 
doi:10. 1088/0954-3899/36/8/085006. 

[30] B. L. G. Bakker, L. A. Kondratyuk, M. V. Terentev, On the for- 
mulation of two-body and three-body relativistic equations employing 
light front dynamics, Nucl. Phys. B158 (1979) 497, doi:10. 1016/0550- 
3213(79)90179-2. 



23 



